Rayleigh- Taylor instability in binary condensates 
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We propose a scheme to initiate and examine Rayleigh- Taylor instability in the two species Bose- 
Einstein condensates. We identify 85 Rb- 87 Rb mixture as an excellent candidate to observe it ex- 
perimentally. The instability is initiated by tuning the 85 Rb- 85 Rb interaction through magnetic 
Feshbach resonance. We show that the observable signature of the instability is the damping of 
the radial oscillation. This would perhaps be one of the best controlled experiments on Rayleigh- 
Taylor instability. We also propose a semi analytic scheme to determinate stationary state of binary 
condensates with the Thomas-Fermi approximation for the axis symmetric traps. 
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Introduction. — Rayleigh- Taylor instability (RTI) sets 
in when lighter fluid supports a heavier one. It is present 
across a wide spectrum of phenomena related to inter- 
face of two fluids. The turbulent mixing in astrophysics, 
inertial confinement fusion and geophysics originate from 
RTI. In superfluids, RTI sets up crystallization waves at 
the superfluid-solid 4 He interface PJ ■ Despite the ubiqui- 
tous nature and importance, controlled experiments with 
RTI are difficult and rare. However, we show that the 
two species Bose-Einstein condensates (TBECs) or bi- 
nary condensates are ideal systems for a controlled study 
of RTI in superfluids. The remarkable feature of TBECs, 
absent in single component BECs, is the phenomenon of 
phase separation. The TBECs, first realized in a mixture 
of two hyperfrne states of 87 Rb are rich systems to 
explore nonlinear phenomena. Several theoretical works 
have examined various aspects of TBECs. These include 
stationary states [!, 0, [f| 0|, modulational instability 
BBIl, collective excitations [H El El El and do- 
main walls solitons[I3|. Another instability related to 
RTI, which has attracted growing interest, is the Kelvin 
Helmholtz instability (KHI). The prerequisites of KHI 
are, phase separation and relative tangential velocities 
at the interface. Quantum KHI have been observed in 
experiments with 3 He [l5[ and recently studied theoreti- 
cally for TBEC [HJ]. 

To initiate RTI we start with the phase separated state. 
Then, increase the scattering length of the species at the 
core. At a certain value it creates a quantum analogue 
of RTI in fluid dynamics. As a case study we choose the 
TBEC of 85 Rb- 87 Rb mixture. In this system, the 85 Rb 
intra species interaction is tunable through a Feshbach 
resonance [l8| and was recently used to study the misci- 
bility [l9j . More recently, the dynamical pattern forma- 
tion during the growth of this TBEC was theoretically 
investigated @. The other feature is, the inter species 
85 Rb- 87 Rb interaction is also tunable and well studied 
[20I ] . Considering the parameters of the experimental re- 
alization, we choose the axis symmetric ( cigar shaped) 
trap geometry. 

Phase separated cigar shaped TBECs. — In the mean 



field approximation, the TBEC is described by a set of 
coupled Gross-Pitaevskii equations 
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where i = 1,2 is the species index, Uu — Aith 2 ai/mi 
with rrii as mass and cij as s-wave scattering length, is 
the intra-species interaction; EZy = 2irK aij /rriij with 
rriij = miTnj/(rn,i + rrij) as reduced mass and as 
inter-species scattering length, is inter-species interac- 
tion and ^ is the chemical potential of the z th species. 
To study the RTI we consider the phase separated state 
(U12 > VU11U22) in axis symmetric trapping potentials 
Vi(p, z) = m i u> 2 (a 2 p 2 + X 2 z 2 )/2. In the present work we 
consider cigar shaped potentials, that is the anisotropy 
parameters oti > Xi and Uij are all positive. Neglecting 
the inter species overlap, the Thomas-Fermi (TF) solu- 
tions are \tpi(p, z)\ 2 = [pi — Vi(p, z)]/Uu. The chemical 
potentials pi are fixed through the normalization condi- 
tions. When oil 3> Aj, the interface of the phase sepa- 
rated state is planar and species having larger scattering 
length sandwiches the other one [l7j |. 

For simplicity of analysis consider trapping potentials 
with coincident centers. Then, let z — ±L\ be the planes 
separating the two components and ±1/2, the spatial ex- 
tent of the outer species along z-axis. The density distri- 
butions n\ and n 2 of the TBEC are 

/ \ Mi - Vi(p,z) 
ni{p,z) = , -Li<z<Li, (2) 
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This assumes no overlap between the two species. Then 
the problem to determine the stationary state is equiv- 
alent to calculating L\. Theoretically, L\ can be deter- 
mined by minimising the total energy of the TBEC with 
fixed number of particles of each species. If Ni and pi 
are the number of atoms and radial size of i th species 
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respectively, then 
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From the TF approximation 
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The total energy of the binary condensate is 

E = [ dV[V 1 (p,z)\^ 1 (p,z)\ 2 + V 2 (p,z)\MP^)\ i 
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We minimize E numerically, with Eq.([5]) and ^ as con- 
straints, to obtain the required value of L±. Substituting 
the value of L\ back into Eqs.© and ([6|), one can de- 
termine pi and p2- Thus Eqs.((5][7| uniquely define the 
stationary state of the TBEC. 

As mentioned earlier, we consider the parameters of 
the recent experiment [19| with 85 Rb and 87 Rb as the 
first and second atomic species. The radial trapping fre- 
quencies are identical (a; = 1) and for the axial trap- 
ping frequencies Ai = 0.022 and A 2 = 0.020. The 
scattering lengths are an = 51a , a 22 = 99a and 
fli2 = 0,21 = 214ao, and we take Ni = 50, 000. Then, 
FigfT] shows the variation in E as a function of L%. The 
value of L\ where minimum of E occurs is 32.5a OS c- Here 
the unit of length a osc = h/niitu with uj = 130Hz, 
is the radial trapping frequency. This is in agreement 
with the numerical result 33.8a osc calculated using split- 
step Crank-Nicholson method (imaginary time propaga- 
tion) We refer to this state as phase I, where 85 Rb 
and 87 Rb are at the center and flanks respectively. We 
have also calculated the equations of interface planes for 
trapping potentials whose minima do not coincide. The 
expressions are much more complicated, however the nu- 
merical and semi-analytic results are in agreement. 

Binary condensate evolution. — In the fluid dynamics 
parlance, the gradient of the trapping potentials are the 
equivalent of gravity. If s is the oscillation frequency of 
the interface between the two condensates, one placed 
over the other. Then from Bernoulli 's p rinciple along 
with proper boundary conditions [22L l23j |. we find from 
linear stability analysis 
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FIG. 1: The variation in energy E with L\ in phase separated 
regime. The upward arrow indicates the position of minimum 
E, which is at L\ — 32.5a osc . Inset shows the same plot along 
with the variation of ^1 and p2 with respect to Li, the blue 
and black curves correspond to pi and pi respectively. 



Here k x and k y are wave numbers along x and y coor- 
dinates. The densities n\ and ri2 are at a point (p, L) 
on the interface. For the sake of simplicity, we con- 
sider mi = TO2 = m and Ai = A2 = A while deriv- 
ing the above relation. There is an instability, referred 
to as Rayleigh- Taylor instability, at the interface when 
m < n 2 . From the TF approximation this condition is 
equivalent to an > a 22 (/ii — V)/(p2 — V). Here V is 
the trapping potential of the two species at the interface. 
Normal fluids with RTI, any perturbation at the inter- 
face however small grows exponentially. Then the lighter 
fluid rises to the top as bubbles and heavier fluid sinks 
as finger like extensions till the entire bulk of the lighter 
fluid is on top of the denser one. On the other hand, 
binary condensates in a similar situation evolve in a very 
different way. 

To examine the dynamical evolution of the binary con- 
densate with RTI, we take phase I ( an < 022) as the 
initial state. In this phase, the 87 Rb BEC at the flanks 
is considered as resting over the 85 Rb BEC at the core. 
Then through the 85 Rb- 85 Rb magnetic Feshbach reso- 
nance [IU increase an till an > a 2 2(pi — V)/(p2 — V) 
to set up RTI. However, maintain U\% > \JUwU22 so 
that the TBEC is still immiscible. Let us call this as 
the phase la and it is an unstable state. The station- 
ary state of the new parameters is phase separated and 
similar in structure to the initial state. But with the 
species interchanged. Let us call the stationary state of 
the new parameters as phase II. The binary condensate 
should dynamically evolve from phase la to II. However, 
unlike in normal fluids with RTI, there are no bulk flows 
of either 85 Rb or 87 Rb atoms, to the periphery of the 
trap. Instead the condensates tunnel with modulations. 
This occurs due to the coherence in the quantum liq- 
uids. To examine the evolution, we solve the pair of 
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time-dependent GP equations 
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which describe the TBEC. During the evolution, the den- 
sity profiles is approximated as rii(p,z) = n° q (p, z) + 
6rii(p,z). Here n° q (p, z) and 5rii(p, z) are the equilib- 
rium density and fluctuation arising from the increase in 
dix- Following the hydrodynamic approximations, the 
5rii{p 1 z) or collective modes follow the equations 
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Srii = Vn, • V UijSrij + n t V 2 ^ Uij6m. (10) 
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Consider Srii(p, z,t) = ai(t)p l exp(±il(f>) as the form of 
the solution, where a^(i) subsumes the time dependent 
part of the solution including temporal variation of the 
amplitude and I is an integer. Then as ~S7 2 5rii — and 
for the miscible phase, considered for simplicity of the 
boundary conditions, we get 
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We can also get a similar set of coupled equations for 
the other form of the collective modes Srii(p, z,t) = 
cii{t)zp l ~ 1 exp(±i(Z — 1)0). In this case the prefactor is 
(I — 1 + X 2 ) instead of I. In either of the cases, the equa- 
tions are similar to two coupled oscillators. For the phase 
separated state, the form of the TF solutions are signif- 
icantly different from the miscible one. However, when 
RTI sets in, the collective modes like in miscible case, 
are damped and coupled as the condensates interpene- 
trate each other. 

TBEC evolution with RTI. — To examine the evolution 
of TBEC with RTI, as mentioned earlier, we choose the 
phase I as the initial state. Then change an to 80ao, 
102ao, 200ao, 306a , 408ao and 780ao, the last value is in 
the miscible parameter region. The dynamical variables 
which are coarse grained representative of the dynamical 
evolution are Prms and z rms , the tuis radial and axial 
sizes. 

When an is increased to 80ao, the 85 Rb condensate os- 
cillates radially to accommodate excess repulsion energy. 
This is the only available degree of freedom as tight con- 
finement, arising from 87 Rb at the flanks, along z-axis 
restricts axial oscillations. In TF approximation the ef- 
fective potential V c g = V+ (p2 — V)U\ilUi2- The angu- 
lar frequency of the oscillation is ss 0.32ui. This is close 
to one of the eigen modes of the Bogoliubov equations. 
The temporal variation of r rms is shown in FigJ5] (inset 
plot). The plots show that, the oscillation of the 87 Rb 
is sympathetically initiated. This is due to the coupling 
between the two condensate species. The oscillations are 
more prominent with less number of atoms. 

There is a change in the nature of oscillations when 
an > 022 (a*i — V)(fjt2 — V). The corresponding station- 
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FIG. 2: The variation in r rms ( in units of a OBC ) for S5 Rb and 
87 Rb with time ( in units of lo" 1 ) when ai is changed from 
51ao to 408ao. The blue and red curves correspond to 85 Rb 
and 87 Rb respectively. 



ary state has 87 Rb and 85 Rb at the core and flank re- 
spectively. The r rms oscillation frequency is the same as 
in an < a 22 case. But there is a temporal decay of the 
amplitude till it equillibrates. The decay is due to the 
expansion of 85 Rb along z-axis and is an unambiguous 
signature of RTI. The expansion is clearly discernible in 
the density profile as shown in Fig[3]and the rate of decay 
increases with Aan. The main plot in Figj2] shows tem- 
poral variation of r Tms for an = 408ao, close to the mis- 
cible domain. There is a strong correlation between the 
decay rate and nature of oscillation. For an marginally 
larger than a22, the 85 Rb condensate tunnels through the 
87 Rb condensate. Where as at larger values the 85 Rb ex- 
pands and spreads into the 87 Rb. 




FIG. 3: Evolution of the TBEC with RTI. The first and sec- 
ond row are density profiles of 85 Rb and 87 Rb BECs respec- 
tively after increasing an to 408ao. Starting from left, the 
density profiles are at 0, 24.5, 49.0 and 73.5 msecs after the 
increase of an. 

A dramatic change of the coupled oscillations occurs 
when U12 < VU11U22, the TBEC is then miscible. 
The 85 Rb expands through the 87 Rb cloud and the two 
species undergo radial oscillations which has a beat pat- 
tern. The FigU shows the r rms when an = 780ao. 
Besides the radial oscillations, as to be expected when 
an > a.22 (Ml — ^)/(a*2 — V) , z lms increases steadily. This 
accommodates the excess repulsion energy along the axial 
direction. Along with the oscillations there are higher fre- 
quency density fluctuations reminiscent of modulational 
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instability. It is to be mentioned that, in earlier works 
0, [H modulational instability in the miscibility domain 
was analysed in depth. For the present case the detailed 
analysis of modulational instability shall be the subject 
of a future publication. 




Summary and outlook. — We have examined the on- 
set of Rayleigh- Taylor instability in TBEC and identi- 
fied the observable signature in the dynamics. We have 
specifically chosen the experimentally well studied 85 Rb- 
87 Rb mixture as case study and propose observing RTI 
with the 85 Rb- 85 Rb Feshbach resonance. Starting from 
an < ct22 5 RTI sets in when the TBEC is tuned to 
an > »22(mi — ^0/(^2 — V) in the TF approximation. 
Then damping of r rms of 85 Rb, species at the core, oscil- 
lations marks the onset of RTI. To analyse the stationary 
states we have proposed a semi analytic scheme, applica- 
ble when A <C 1, to minimize the energy functional with 
TF approximation. The results of which are in excellent 
agreement with the numerical results. The A <C 1 is also 
the case when the interface is planar and RTI is more 
prominent. 
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FIG. 4: The variation in r rms ( in units of a OBC ) for 85 Rb and 
87 Rb with time ( in units of when ai is changed from 

51ao to 780ao. The blue and red curves correspond to 85 Rb 
and 87 Rb respectively. 
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